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Abstract 

We study dilepton production from pp, pd and pBe collisions from 1-5 GeV in- 
cluding the 7T°, 77, ijj and A Dalitz decays, direct decays of vector mesons (p, uj) as well 
as subthreshold p production via baryonic resonances (e.g. Z?i3(1520), P\\ (1710)). 
Our calculations compare rather well with the pp and pd data from the DLS Collab- 
oration, however, overestimate slightly the 'old' pBe data from that group. Futher- 
more, detailed predictions for differential dilepton spectra at SIS energies are made 
with a high mass resolution that can be controlled experimentally by the HADES 
Collaboration in near future. 
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1 Introduction 



Electromagnetic probes such as dileptons provide the most clear signals from the early 
phases of high-energy heavy-ion collisions because they may leave the reaction volume es- 
sentially undistorted by final-state interactions. Differential dilepton spectra from heavy- 
ion collisions thus can provide information about the effective degrees of freedom at high 
baryon density and temperature. Apart from that the in-medium properties of hadrons 
(cf. Refs. |1|, 0, |3|, [|, U), that are closely related to the latter problem, can be ex- 
plored as well. According to QCD sum rules 0, [|, [7| as well as QCD inspired effective 
Lagrangian models JTJ, ^, [|, ^ |], [TO], [□], |T2|, |T3|, [T4| especially the vector mesons (p, uj and 
4>) significantly change their properties with nuclear density. It has been argued, that the 
experimentally observed enhanced production of soft lepton pairs above known sources in 
nucleus-nucleus collisions at SPS energies [T5|, |TB| might be due to the in-medium modifi- 



cation of vector mesons [|17], [Tj| [19], £1]] rather than reflecting a new state of hadronic 
matter. 

Dileptons from heavy-ion collisions have also been measured by the DLS Collabora- 
tion p2l E3| at BEVALAC energies, where a different temperature and density regime is 



probed. The first generation of DLS data |22|], based on a limited data set, were consis- 
tent with the results from early transport model calculations p|, ER EBI, I27|] including the 



conventional dileptons sources as pn bremsstrahlung, n°, r], uj and A Dalitz decay, direct 
decay of vector mesons and pion-pion annihilation, without incorporating any medium 
effects. A more recent DLS measurement P3| including the full data set and an im- 
proved analysis, however, shows an increase by about a factor of 5-7 in the cross section 



in comparison to Ref. |22| and the early theoretical predictions. 

This enhancement of the e + e~ yield in nucleus-nucleus collisions provides some theoret- 
ical 'puzzle'. In Ref. [28| the in-medium p spectral functions from Refs. [|1^, [14] have been 
implemented in the HSD transport approach for p mesons produced in baryon-baryon, 
pion-baryon collisions as well as from irir annihilation, and a factor of 2-3 enhancement 
has been obtained compared to the case of a free p-spectral function. In Ref. drop- 
ping vector meson masses and uj meson broadening were incorporated in the UrQMD 
transport model, which also showed an enhancement of the dilepton yield, however, the 



authors could not describe the new DLS data [23] as well. Another attempt to solve the 
DLS 'puzzle' has been performed in Ref. |3(| where the dropping hadron mass scenario 
was considered together with the subthreshold p production in ttN scattering via the 



baryonic resonance iV"(1520) whose importance was pointed out in Refs. JTJ, |31| . It was 
found that the enhancement of the dilepton spectra due to low mass p's from the iV(1520) 
was not sufficient to match the DLS data. Thus, all in-medium scenarios that successfully 
have explained the dilepton enhancement at SPS energies failed to describe the new DLS 
dilepton data ]23| from heavy-ion collisions ( 12 C + 12 C and 40 Ca + 40 Ca). 
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In 1998 the DLS Collaboration has, furthermore, published dilepton data from elemen- 



tary pp and pd collisions at 1-5 GeV |32| which provides the possibility for an independent 



check of the elementary dilepton channels that enter as 'input' in transport calculations 



for heavy-ion reactions. Such an analysis has been carried out in Refs. [££], |33| and it was 
shown that the dilepton invariant mass spectra from pp reactions in the energy range from 
1-5 GeV could be rather well described within the standard sources of hadronic Dalitz 
and direct decays when including also the 'subthreshold' production of p-mesons by bary- 
onic resonances Q. It has been, furthermore, argued that additional channels (like dilep- 
ton decays of scalar mesons (ao, fo) or four-body final states decays like rj' — > 7r + 7r~e + e~ , 
(p°,uj) — > 7r + 7r~e + e - , (p°,uj) — > 7]7T°e + e~) might contribute to the elementary yield in pp 
and pd reactions f34j. Thus it is indispensible, that also these elementary reactions have 
to be measured again with high resolution in order to understand more accurately the 
e + e~ radiation mechanisms. 

In this article we perform a detailed study of dilepton production from pp, pd and 
pBe collisions from 1-5 GeV including all channels as described in Ref. fl3"3"| , i.e. the 
subthreshold p production via baryonic resonances (e.g. D 13 (1520), P n (1710)) as well 
as 7T°, r], lu and A Dalitz decays and the direct decay of vector mesons (p, u, <fi). We 
investigate the role of isospin degrees of freedom with respect to dilepton production in 
pd reactions by employing exact energy-momentum conservation for the deuteron target. 
The pBe reactions are calculated within the resonance transport model [^] that includes 
the same elementary production amplitudes. 

The paper is organized as follows: In Section 2 we briefly present the underlying res- 
onance model that enters the coupled-channel BUU transport approach. In Section 3 we 
provide a comparison of our results with the corresponding DLS data including the DLS 
acceptance and resolution functions. In Section 4 we, furthermore, give detailed predic- 
tions for these reactions with a high mass resolution in view of upcoming experiments 
with the HADES detector at GSI Darmstadt. We close with a summary and discussion 
of open problems in Section 5. 

2 Description of the model 

2.1 Resonance approach 

We perform our analysis of dilepton production from pp, pd and pBe collisions within 
the resonance approach of Refs. [135], |3"6|. This model is based on the resonance concept 



of nucleon-nucleon and meson-nucleon interactions at low invariant energy ^fs |37j by 



adopting all resonance parameters from the Manley analysis [38]; all states with at least 



1 For a detailed description of the individual channels we refer the reader to Ref. 
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2 stars in Ref. gS] are taken into account: P 33 (1232), P n (1440), D 13 (1520), 5 U (1535), 
P 33 (1600), S 31 (1620), 511(1650), £>i 5 (1675), F 15 (1680), Pi 3 (1879), S 31 (1900), P 35 (1905), 
P 3 i(1910), P> 35 (1930), F 37 (1950), Fi 7 (1990), G 17 (2190), P> 35 (2350). These resonances cou- 
ple to the following channels: iVvr, Nrj, Nuj, AK, A(1232)tt, Np, Na, iV(1440)7r, A(1232)p 
with respect to the production and decay. 

It has been shown, that the resonance model provides a good description of the exper- 
imental data on one- and two-pion production in nucleon-nucleon collisions at low energy 
J37|. However, with increasing energy the resonance contributions underestimate the 



data; the missing yield is then treated as a backgraund term to the resonance amplitude. 
This background term 'mimics' t-channel particle production mechanism as well as other 
non-resonance contributions (e.g., direct NN — > NNtt, without creating an intermediate 
resonance). 

With increasing energy, furthermore, the multiparticle production becomes more and 
more important. The high energy collisions - above y/s = 2.6 GeV for baryon-baryon 
collisions and \fs = 2.2 GeV for meson-baryon collisions - are described by the LUND 
string fragmentation model FRITIOF [p9| |. This aspect is similar to that used in the HSD 
approach ||, H, gl], g2| and the UrQMD model @. 

This combined resonance-string approach allows to calculate particle production in 
baryon-baryon and meson-baryon collisions from low to high energies. The collisional 
dynamics for proton-nucleus reactions, furthermore, is described by the coupled-channel 
BUU transport approach [BR FJH] that is based on the same elementary cross sections. 



2.2 Dilepton production 

The dilepton production within the resonance model can be schematically presented in 
the following way: 

(1) 
(2) 
(3) 

e + e~X, (4) 
e + e~X. (5) 

In words: in a first step a resonance R might be produced in baryon-baryon (BB) or 
meson-baryon (mB) collisions - (0), (|2]). Then this resonance can couple to dileptons 
directly - @ (e.g., Dalitz decay of the A resonance: A — > e + e~iV) or decays to a meson 
m (+ baryon) - (f|) which produces dileptons via direct decays (p, to) or Dalitz decays 
(tt° , 7] , uj) . The resonance R might also decay into another resonance R' - (^) which 
later produces dileptons via Dalitz decay or again via meson decays (e.g., P 35 (1930) — > 



BB RX 

mB RX 

R e + e-X, 
R — > mX, m - 
R — > R ' X } R' 
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Ap, A — > e + e _ iv, p — > e + e~). Note, that in our combined model the final particles 
- which couple to dileptons - can be produced also via non-resonant mechanisms, i.e. 
'background' at low and intermediate energies and string decay at high energies. 

The electromagnetic part of all conventional dilepton sources - 7T°,r],uj and A Dalitz 
decay, direct decay of vector mesons p,ui - are treated in the same way as described 
in detail in Sec. 3 of Ref. |33[] where we have calculated the dilepton spectra from pp 
collisions in comparison to the DLS data. The only difference is that in our present 
calculation we use the spectral function of vector mesons (p, u) in the same form as for 
all other resonances as used in Manley analysis ||38|| : 



2 M 2 Y V 

f(M) = N V - — v , 2 - (6) 

7T (M 2 - mi,) 2 + {MTlt} 2 

In (H) Ny guarantees normalization to unity, i.e. / f(M)dM = 1, while TY ot is the total 
vector meson width. 

For pd and pBe collisions we also take into account pn bremstrahlung using the soft- 
photon approximation ||44|| , where the radiation from internal meson lines is neglected and 
the strong interaction vertex is assumed to be on-shell. In this case the strong interaction 
part and the electromagnetic part separate; however, the cross section for dileptons has to 
be corrected P5| , (4"6|] by reducing the phase-space for the colliding particles in their final- 
state. We discard pp bremsstrahlung from the nucleon pole since the microscopic OBE 
calculations j|7|, ^8], ^9f have shown that at 1.0 GeV the pp bremsstrahlung is smaller than 
the A Dalitz decay contribution by a factor of 2-3. At this energy also the interference 
between nucleon and A-pole terms is negligible. At high energies, where this interference 
becomes important ||7]], however, the overall contribution from these channels is negligible 
(see below). We, furthermore, discard e + e~ production from the scalar mesons fo,a as 
advocated in Ref. [34| since our studies in Ref. |50| have indicated that the multiplicity 
of scalar mesons in pp reactions is very small in the energy regime of interest. 

The basic underlying assumption for dilepton production, that enters in our calcula- 
tion, is that the amplitude can be factorized into a meson or/and resonance production 
and a dilepton decay part. This assumption provides not only a significant simplification, 
which is necessary for applications in transport calculations [25, 42, 35], but also allows to 
take into account the 'inelasticity' from many-particle production channels which become 
dominant at high energy. 

In Fig. [I] we illustrate our factorization assumption for the process NN — > RN — > 
pNN — > e + e~NN with an intermediate resonance R coupled to the p meson, e.g. _Di3(1520), 
which is responsible for p production below the threshold. As shown in Fig. [I] we cut the 
diagram twice such that we have a resonance production part (NN — > RN), a resonance 
decay to the p meson (R — > pN) and direct dilepton decay of the p meson (p — > e + e~). 
This assumption with respect to subthreshold p production via the _Di3(1520) is different 
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to the ansatz used in the previous dilepton analysis [[!(], [33|] . In Ref . we have factor- 
ized the amplitude NN — > RN pNN — > e + e~NN into a _D 13 (1520) production part 
(NN -> RN) and Dalitz decay part R -> piV -> e+e~iV (cf. Sect. 3.1 of Ref. In 
the latter approach the p meson is considered to be off-shell, i.e. a virtual particle, and 
can have a mass below 2m n which leads to the Dalitz behavior of the dilepton spectra 
(cf. Fig. 3 of Ref. |3~3||), whereas in the present approach we assume that the p meson 
has a free spectral function with masses only above 2m n . The latter assumption leads to 
the 'direct p decay' profile for dileptons from subthreshold p's and is consistent with our 
resonance prescription used also in the transport calculations where we propagate 
such p mesons explicitly. 

To compare with the experimental data of the DLS Collaboration one has to use the 
appropriate experimental filter, which is a function of the dilepton invariant mass M, 
transverse momentum p^ and rapidity in the laboratory frame - F (M , p^ , Uub) • F° r 
that purpose we explicitly simulate all dileptons by Monte-Carlo and apply the filter 
F(M,p T , 



3 e + e results in comparison to DLS data 
3.1 Dileptons from pp collisions 

In Fig. |2| we present the calculated dilepton invariant mass spectra da/dM for pp collisions 
from 1.0 - 4.9 GeV including the final mass resolution and filter F(M,pT,yi a b) from the 
DLS Collaboration (version 4.1) in comparison to the DLS data f32fl . The thin lines 
indicate the individual contributions from the different production channels; i.e. starting 
from low M: Dalitz decay 7r° — > 7e + e~ (short dashed line), 77 — > 7e + e~ (dotted line), 
A — > Ne + e~ (dashed line), uj — > 7r°e + e~ (dot-dashed line), for M w 0.7 GeV: to — > e + e~ 
(dot-dashed line), p° — > e + e~ (short dashed line). The full solid line represents the sum 
of all sources considered here. 

Whereas at 1.04 GeV the dileptons stem practically all from 7r° and A Dalitz decays, p 
decay becomes more important with increasing energy and at high invariant mass. Note, 
that in the present analysis we do not separate the production mechanism for p mesons 
as in Ref. p3| , i.e. the p- meson contribution might come from the decay of Z?i3(1520) 



and Pn(1710) or non-resonance p production. Thus, the short dashed line (denoted as 
p° — > e + e~) presents the sum of all p contributions. However, the resonance mechanism 
dominates for p production at low energies, especially below the pp — > ppp threshold 
y/s < 2rri]\T + m p , where yfs is the invariant collision energy (see Sec. 2 in Ref. f33j). 



At 2.1 and 4.9 GeV the dilepton yield is dominated by the 77 Dalitz decay and direct 
decays of p and uo mesons. The p spectrum is enhanced towards low M due to the limited 
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phase space and strong mass dependence (M~ 3 ) of the dilepton decay width |33] in line 
with the vector dominance model. Note, that at 2.1 GeV we use the resonance production 
mechanism in our model, whereas at 4.9 GeV most of the particles are produced by string 
decays. Within the strict resonance model which gives only the exclusive pp — > ppr] cross 
section the r\ yield from pp collisions is underestimated by about a factor of 2 (dot-dot- 
dashed line) at 2.1 GeV. Since at this energy other channels such as NNrjn, NNrjixir are 
open, we have estimated the cross section for these channels by the string model and 
added it to the total r\ yield. This gives the final rj contribution (dotted line) for 2.1 GeV. 



3.2 Dileptons from pd collisions 

The experimental measurement of dilepton spectra from proton and deuteron targets 
provides a unique possibility to get information about the isospin dependence of the 
NN interaction which is interesting in itself and, furthermore, important with respect to 



heavy-ion reactions. The pp and pd dilepton data from the DLS Collaboration [32] are 
the first experimental step in this direction. 

Since the deuteron is a bound system of proton and neutron, the Fermi motion of 
the nucleons plays an important role for all interactions. For example, Fermi motion 
leads to a shift of the threshold for particle production in pd interactions to lower y/s, 
i.e. the channels closed by kinematics in pp collisions might be open in pd. Thus, the 
pd reactions can not be directly approximated by the simple sum of pp and pn. On the 
other hand, one can't treat the deuteron within the usual mean-field approach used in 
transport calculations for more heavy nuclei. 

In our analysis we thus take into account the Fermi motion of nucleons using the 
momentum distribution of the deuteron calculated with the Paris potential ||51|| . We 



initialize a deuteron in momentum space (using a Monte Carlo method) according to the 
momentum distribution f(p): 

m = (7) 



where p is the momentum of a nucleon in the deuteron center-of-mass, while ^(p) is the 
deuteron wave function. To fulfill the binding energy constraint we assume 



E spec — y p 2 + tti 2 N} E part — E$ — E spec . (8) 

The nucleon, which does not participate in the reaction, i.e. the spectator, has a free 
dispersion relation for the energy-momentum E spec , whereas E part is the energy of the 
nucleon participating in the reaction (participant). In (H) Ed = 2mjv + e is the total 
deuteron energy, where e = —2.2 MeV is the binding energy of the deuteron. 

The invariant collision-energy distribution dN/dyfs for pd reactions calculated with 
the dispersion relation (|8]) is shown by the solid histogram in Fig. |3|. For comparison 
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(dashed histogram) we present the dN/dy/s distribution in the physical region (2m N < 
\fs < y/s max) calculated with f(p) (0) in the quasi-free scattering limit, i.e. using the free 



dispersion relation for both nucleons in the deuteron: E spec = E part = \Jp 2 + mj v . The 
lower limit of y/s is given by the free cross section that only exists for \Z>2m^. As seen 
from Fig. |||, both prescriptions are the same in the peak value, however, the deuteron 
spectral distribution (^) leads to an enhancement of the y/s distribution at small y/s and 
a reduction at high y/s. This implies a substantial lowering of particle production at high 
y/s with respect to the quasi- free limit, which has been often used in the literature. 

Within the 'binding' prescription ([8]) we have calculated the dilepton invariant mass 
spectra da/dM for pd collisions from 1.0 - 4.9 GeV including the final mass resolution 
and filter F(M,pT,yi a b) from the DLS Collaboration (version 4.1). The results are show 
in Fig. f| in comparison to the DLS data P^| , where the assignment of the lines is the 
same as in Fig. |2|. Additionally, the dot-dot-dashed lines indicate the pn bremsstrahlung 
contributions. 

Similar to Fig. |2|, at 1.04 GeV the ir° and A dalitz decays are dominant, however, 
due to the larger mass of the deuteron, the n channel is already open contrary to the 
pp reaction and the subthreshold p branch is getting more important at high invariant 
mass M. Also pn bremsstrahlung contributes at low energies, however, with increasing 
bombarding energy it becomes less important. 

Starting from 1.61 GeV the dominant contribution at low M is again the 77 Dalitz 
decay, whereas at high M the dileptons stem basically from direct p, lu decays. The 77 
contribution here is calculated by assuming a pn ^ pnr] = 6cr pp _> ppri for 2.425 GeV < y/s < 



2.517 GeV in line with the experimental data from Ref. ||52|| , o pn ^ pnr] = 2.5a pp ^ ppr] for 
2.517 GeV < y/s < 2.65 GeV and <J pn ^ pnr] ~ cr pp ^ pprj for y/s > 2.65 GeV where the 
multiparticle production (above 2.65 GeV) is evaluated within the LUND string model. 
Note, that at high y/s the isospin asymmetry a pn /a pp should be washed out by the isospin 
of pions that are created together with the rj. Though this 'recipe' seems to work in 
comparison to the DLS data, we point out that further experimental data with good 
mass resolution will be necessery to clarify the situation for the isospin couplings. 

In order to demonstrate the importance of isospin degrees-of-freedom for dilepton 
production we present in Fig. [| the invariant mass spectra da/dM for pp (short dashed 
line), pn (long dashed line), pp + pn (dot-dashed line) and pd (solid line) collisions at 1.61 
GeV in comparison to the DLS data [5^] . At these energies the pn channel gives the main 
contribution to the pd dilepton yield. Furthermore, the influence of the Fermi motion 
in the deuteron is clearly seen at high invariant mass M by comparing the dot-dashed 
(pp + pn) and solid (pd) lines. 

In Fig. P we present the excitation function, i.e. the integrated dilepton cross section 
for masses above 0.15 GeV, as a function of the available energy Q = y/s — 2m p for pp 
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(dashed line) and pd (solid line) in comparison to the DLS data (open circles for pp 
and full circles for pd). The calculated excitation functions for pp and pd increase with 
the energy and stay within the experimental errorbars. 



3.3 Dileptons from pBe collisions 

Since Be is one of the lightest nuclei, in-medium effects play a minor role because most 
hadron decays occur already in vacuum. So one can consider pBe collisions as an addi- 
tional and independent check for the elementary reactions and isospin degrees of freedom. 
However, the collisional dynamics has to be taken into account properly. For that purpose 
we use the coupled-channel BUU transport approach as discussed in Sec. 2. 

Dileptons from pBe have been measured by the DLS Collaboration in the end of the 
80's f2"2| . In Fig. we show the calculated dilepton invariant mass spectra da/dM for 
pBe collisions at 1.0 (upper part), 2.1 (middle part) and 4.9 GeV (lower part) including 
the final mass resolution and filter F(M,pT,yi a b) from the DLS Collaboration (version 



1.6) in comparison to the data [22]. The assignment of the lines is the same as in Figs. |^ 
and §. With the filter 1.6 we do not get enough suppression for the dilepton yield at 
M < 0.15 GeV which basically stems from tt° Dalitz decay. Similar to pd reactions the 
dominant channels at 1.0 GeV are the A Dalitz decay and pn bremsstrahlung. At 2.1 and 
4.9 GeV the 77 Dalitz decay becomes important at low M and the direct decays of vector 
mesons (p, uj) for M ~ 0.8 GeV. Within the present model we slightly overestimate the 
data at 4.9 GeV, however, we discard any discussion about the physical reason since the 
'old' DLS filter is not sufficiently understood. For example, the 'old' DLS data for pBe at 
4.9 GeV show a similar dilepton yield as for pd, whereas according to our calculations it 
should be higher by a factor of ~ 2 since more nucleons are involved in the pBe reaction. 

In this respect we have performed detailed calculations for pp, pd and pBe reactions 
without implementing any filter, however, including a high (10 MeV) mass resolution for 
da/dM. For the results, that can be controlled by the HADES Collaboration in near 
future, we refer to the following Section. 



4 Predictions for HADES 

In this Section we present our predictions for the dilepton invariant mass spectra, trans- 
verse momentum (p?) distributions and rapidity (y) distributions for pp, pd and pBe 
collisons at 1.5 GeV and 4.0 GeV, which will be measured with the HADES spectrometer 
at GSI. 

In Fig. [| the calculated dilepton invariant mass spectra da/dM for pp (upper part), 
pd (middle part) and pBe collisions (lower part) at 1.5 GeV (left panel) and 4.0 GeV 
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(right panel) are shown including a AM = 10 MeV mass resolution. The assignment of 
the lines is the same as in Figs. |2| and [|. The dominant contribution at low M (> m n o) 
is the 77 Dalitz decay, however, for M > 0.4 GeV the dileptons stem basically all from 
direct vector meson decays (p and u). Within a 10 MeV mass resolution the peak from 
the direct uj decay (which was completely smeared out within the DLS mass resolution) 
is seen. The u peak is getting more pronounced at high energy (cf. the right panel of 
Fig. |). 

In Fig. H we present the transverse momentum distribution da / dp? / (2tcpt) for pp 
(upper part), pd (middle part) and pBe collisions (lower part) at 1.5 GeV (left panel) and 
4.0 GeV (right panel). The assignment of the lines again is the same as in Figs. ^| and 
[|. Here we did not imply any cut in the invariant mass M. For all systems and energies 
considered the tt° Dalitz decay is much higher than all other contributions, which is in 
line with the result of Fig. [8] because the pion production cross section is much larger 
than that for the other mesons. The shape of the ir°,p,A,uj Dalitz decays is defined by 
the kinematics of the Dalitz decay and the 'flat' part at low pt depends on the mass of 
the decaying resonance. 

The laboratory rapidity distributions da/dy for pp (upper part), pd (middle part) and 
pBe collisions (lower part) at 1.5 GeV (left panel) and 4.0 GeV (right panel) are shown in 
Fig. pi], where the assignment of the lines is the same as before. The situation is similar 
to Fig. - the 7T° Dalitz decay is the only visible contribution when imposing no cuts on 
the invariant mass. 

In order to make the contribution from other channels visible one has to suppress the 
'background' tt° yield. In Fig. |Tl| we display the transverse momentum spectra for pp 
(upper part), pd (middle part) and pBe collisions (lower part) at 1.5 GeV (left panel) 
and 4.0 GeV (right panel) using a cut in invariant mass as 0.4 GeV< M < 0.7 GeV. 
Fig. [12| shows the laboratory rapidity distributions for the same systems and kinematical 
conditions. With this cut on intermediate dilepton masses one can select the region where 
the 7T° Dalitz decay is totally removed and the r\ Dalitz decay is suppressed substantially 
(cf. Fig. H); in this way the p branch in the transverse momentum and rapidity spectra 
becomes seen (cf. Figs. [IT], |TT|). Thus, a measurement of the transverse momentum and 
rapidity spectra with a corresponging cut in M will provide additional and independent 
information about the p contribution to the dilepton spectra. 

Finally, in Fig. [13| we present the double differential dilepton spectra da/dMdp? as 
a function of the invariant mass M and transverse momentum p? for pd collisions at 4.0 
GeV. At fixed pt one can recognize the shape of the invariant mass spectra (cf. Fig. H) 
with a strong ir° Dalitz decay branch at low M as well as the contributions from 77 Dalitz 
decay and the vector meson (p, ui) decays. At fixed M the shape looks similar to the 
one in Fig. [||. At low M the exponential decrease stems from tt° Dalitz decay, then the 
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spectra become flatter due to the contributions from 77 Dalitz (M < 0.4 GeV) and direct 
p decays. At M ~ 0.78 GeV the peak from the direct decay of u mesons is visible. Thus, 
such type of 3-dimensional experimental information (or even 4-dimensional including 
rapidity) allows to select the contributions from different channels. 

5 Summary 

Within the framework of the combined resonance-string approach we have performed a 
detailed study of dilepton production from pp, pd and pBe collisions at 1.0 - 4.9 GeV. In 
this approach the particle production proceeds through resonances at low energies and by 
string formation and decay at high energies. The collisional dynamics is described by the 
coupled-channel BUU transport approach. As dilepton sources we have considered n°, r], 
uj and A Dalitz decays, pn bremsstrahlung and direct decays of vector mesons (p, uj). 

We have compared our calculated results to the experimental data of the DLS Col- 
laboration which provide constraints on the different individual production channels. It 
has been found that the DLS data for the pp and pd collisions can be reasonably well 
described in our model whereas for pBe systems (especially at 4.9 GeV) our calculations 
give a slightly higher dilepton yield. We have demonstrated the importance to measure 
the dileptons from pp and pd (or even pBe) collisions simultaneously since such data 
can provide constraints on the isospin dependence of pp and pn interactions, which is 
important for an understanding of heavy-ion data. 

Dilepton spectra from elementary reactions will be measured in future by the HADES 
detector at GSI Darmshtadt with high mass resolution and good accuracy. In this respect 
we have made predictions for the dilepton invariant mass spectra, transverse momentum 
and rapidity distributions for pp, pd and pBe at 1.5 and 4.0 GeV. We showed that proper 
cuts in invariant mass for transverse momentum and rapidity spectra allow to select 
different dilepton sources and to study, for example, the p meson channel in more detail. 
We have indicated, furthermore, that it might be very useful to provide experimentally 
multi-dimensional information, e.g. double differential dilepton spectra da/dMdpT, in 
order to investigate the individual contributions. 
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15 




Figure 2: The calculated dilepton invariant mass spectra da/dM for pp collisions from 
1.0 - 4.9 GeV in comparison to the DLS data [32]. The thin lines indicate the individual 
contributions from the different production channels; i.e. starting from low M: Dalitz 
decay n° — > •ye + e~ (short dashed line), rj — > •ye + e~ (dotted line), A — ► Ne + e~ (dashed 
line), u — > 7r°e + e~ (dot-dashed line); for M ~ 0.7 GeV: uj —> e + e~ (dot-dashed line), 
p° — > e + e~ (short dashed line). The full solid line represents the sum of all sources 
considered here. At 2.09 GeV the dot-dot-dashed line indicates the r\ contribution from 
the resonance decays only while the dotted line corresponds to the inclusive r\ production 
(see text). 
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Figure 3: The distribution in the invariant collision energy y/s {dN/d\fs) in arbitrary 
units for the pd reaction at 1.04 GeV calculated with the dispersion relation (||) (solid 
histogram) and within the quasi- free scattering limit (dashed histogram). 
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Figure 4: The dilepton invariant mass spectra da/dM for pd collisions from 1.0-4.9 GeV 
in comparison to the DLS data |32|. The assignment of the lines is the same as in Fig. 
Additionally, the dot-dot-dashed lines indicate the pn bremstrahlung contributions. 
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M [GeV] 

Figure 5: The dilepton invariant mass spectra da/dM for pp (short dashed line), pn 
(long dashed line), pp + pn (dot-dashed line) and pd (solid line) collisions at 1.61 GeV 
in comparison to the DLS data |32fl . The difference between the dot-dashed line and the 
solid line stems from the Fermi motion of the nucleons in the deuteron. 
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Figure 6: The excitation function for dileptons from pp (dashed line) and pd (solid line) 
in comparison to the DLS data |32| (open circles for pp and full circles for pd reactions). 



20 





M [GeV/c 2 ] 



Figure 7: The dilepton invariant mass spectra da/dM for pBe collisions at 1.0 GeV (upper 
part), 2.1 GeV (middle part) and 4.9 GeV (lower part) in comparison to the DLS data 
p2| . The assignment of the individual lines is the same as in Figs. § |3| 
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Figure 8: The dilepton invariant mass spectra da/dM for pp (upper part), pd (middle 
part) and pBe collisions (lower part) at 1.5 GeV (left panel) and 4.0 GeV (right panel) 
including a 10 MeV mass resolution. The assignment of the individual lines is the same 
as in Figs. ||. f|. 
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p T [GeV/c] p T [GeV/c] 

Figure 9: The transverse momentum distribution da / dp? / (2itpt) for pp (upper part), pd 
(middle part) and pBe collisions (lower part) at 1.5 GeV (left panel) and 4.0 GeV (right 
panel). The assignment of the individual lines is the same as in Figs. |j. 
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Figure 10: The rapidity distributions da/ dy for pp (upper part), pd (middle part) and pBe 
collisions (lower part) at 1.5 GeV (left panel) and 4.0 GeV (right panel). The assignment 
of the individual lines is the same as in Figs. 0, [|. 
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Figure 11: The transverse momentum distribution da/dpT/ (2ttpt) for pp (upper part), pd 
(middle part) and pBe collisions (lower part) at 1.5 GeV (left panel) and 4.0 GeV (right 
panel) implying a cut in invariant mass of 0.4 < M < 0.7 GeV. The assignment of the 
individual lines is the same as in Figs. 0, £1 



25 



0.4 < M < 0.7 GeV 




Figure 12: The laboratory rapidity distributions da/dy for pp (upper part), pd (middle 
part) and pBe collisions (lower part) at 1.5 GeV (left panel) and 4.0 GeV (right panel) 
implying a cut in invariant mass of 0.4 < M < 0.7 GeV. The assignment of the individual 
lines is the same as in Figs. |2|, |j. 
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Figure 13: The calculated double differential dilepton spectra da / dMdpx as a function of 
invariant mass M and transverse momentum pt for pd collisions at 4.0 GeV. 
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